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Abstract 

We will give an outline of the computation of the QCD corrections to the spin structure 
function gi(x,Q 2 ) and the spin fragmentation function g^(x,Q 2 ) which are measured in 
deep inelastic electron-proton scattering and in electron-positron annihilation respectively. 
In particular we show how to deal with the 75-matrix and the Levi-Civita tensor, ap- 
pearing in the amplitudes of the parton subprocesses, when the method of iV-dimensional 
regularization is used. 
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Figure 1: Kinematics of polarized deep inelastic electron-proton scattering 



1 Deep inelastic electron-proton scattering 

Deep inelastic electron-proton scattering proceeds via the following reaction (see Fig. [1]) 

e-(h, aO + P(p, s) -> e-(l 2 , cr 2 ) +' X' . (1) 

Here 'X' denotes any inclusive hadronic final state and V in Fig. [I] stands for the neutral 
intermediate vector bosons given by j,Z. For simplicity we will assume that the momentum 
transfer is very small with respect to the mass of the Z-boson so that the process in Fig. [I] is 
dominated by the one photon exchange mechanism only. In the case the proton is polarized 
parallel (— >) or anti-parallel (<— ) with respect to the spin of the incoming electron we obtain the 
cross section 

d 2 a(^) d 2 a(^) Ana 2 r 
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where gi(x, Q 2 ) denotes the longitudinal spin structure function. Further we have defined the 
scaling variables 

Q 2 



X 



y 



p-q 



Q 



-Q < . 



2p ■ q p ■ h 

The spin structure functions show up in the antisymmetric part of the hadronic tensor 
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which is given by 



m 
2p ■ q 



p-q 



(3) 



(4) 



(5) 



Here g2(x,Q 2 ) denotes the transverse spin structure function which is kinematically suppressed 
in cross section (|2]). Since the leading power corrections are of twist two, one can give a parton 
model description of the longitudinal structure function which can be written as 
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Unfortunately there does not exist such a simple formula for g2(x,Q 2 ) because of twist three 
contibutions which are not power suppressed with respect to the twist two parts. Hence one 
cannot give a simple parton model interpretation for the transverse spin structure function and 
we will therefore not discuss it in the subsequent part of this paper. 

In Eq. (§) we have used the following notation. The charge of the light quarks is denoted by e k 
and rif stands for the number of light flavours. The spin parton densities Afi(z,[i 2 ) (i = q,g) 
depend on factorization scale \x which is put to be equal to the renormalization scale. The spin 
parton coefficient functions AC\^ depend on the same scale \i. The quark parton densities and 
the quark coefficient functions can be split in non-singlet (NS) and singlet (S) parts with respect 
to the flavour group. The singlet and non-singlet combinations of parton densities are given by 

"/ 

A/ ? W 2 ) = £[A/ fc (^/i 2 ) + A/ s (z,/i 2 )], (7) 

k=l 

and 

A/, N fc s (z,/i 2 ) = Af k (z^ 2 ) + Af- k (z,fi 2 ) - —Af^z^ 2 ) , (8) 

n f 

respectively. Since it turns out that the equations are easier to study when one performs a Mellin 
transform defined by 

F (n) = ( l dzz n - l F(z) , (9) 



we will present all the following formulae in this representation. 

The parton densities and the coefficient functions above satisfy the renormalization group equa- 
tions. Let us first define the differential operator 

c\ c\ 3 

D = ^— + (3{g)— , 0k) = _#,-!_ + ... . (10) 
d a/i og 1D7T 

Using this notation the renormalization group equations for the parton densities read 

DA/ NS,(n) = _ A7 NS,(n) A/ NS,(n) k = u , d ... 

DAft>M = -At&^AZ/'W , i,j = q,g , (11) 
and for the coefficient functions 

jDAC NS,(n) = A7 NS,(n) AC NS,(n) 

DAC!j a) = A^ACff , i,j = q,g . (12) 

From the equations above it follows that the structure function is a renormalization group 
invariant i.e. 

Dgi n \Q 2 ) = 0, (13) 

which implies that it is a physical quantity independent of the scale \x. 

The anomalous dimensions and the coefficient functions are calculable order by order in 
perturbation theory. Let us first sketch the derivation of the anomalous dimensions before we 
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pay attention to the coefficient functions. In |]J the anomalous dimensions appearing in the spin 
dependent quantities have been derived from the calculation of the operator matrix elements 
(OME's). For an alternative derivation see ||. These OME's are obtained by sandwiching local 
operators between quark and gluon states. These operators appear in the lightcone expansion 
of the product of the electromagnetic currents in Eq. (Q). Suppressing some irrelevant Lorentz 
indices the expansion reads as follows 

oo 

^^rEEfl^C) > i = 9,g • (14) 

Here n denotes the spins of the local operators Of 1 ^ and c[™^(x 2 ) are the Fourier transforms of the 
coefficient functions in position space. The operators of twist two, which can also be split into 
singlet and non-singlet parts, are given in the literature (see e.g. Eqs. (2.5)-(2.7) in [|TJ). Since, 
in the Bj0rken limit, the integrand in Eq. (|J) is dominated by the light cone we can replace the 
current-current product by the above expansion so that one has to compute the renormalized 
OME's 

At\ z ^ r ) = m\ot ) \vm, (15) 

where i, j — q,g and p denotes the external momentum of the quarks and gluons and r = NS, S. 
In and recently also in || the above operator matrix element (OME) has been computed 
up to second order in the strong coupling constant a s . The calculation proceeds as follows. 
After having derived the operator vertices (see Appendix A in Jl]) one has to compute the 
Feynman graphs (see Figs. 1-6 in [l[]) which correspond to the unrenormalized (bare) OME's. 
The latter reveal ultraviolet singularities which are regularized by N-dimensional regularization. 
The unrenormalized OME's indicated by a hat can be written in the form 

^ (n) p£,-) = 5 tl + a s S £ (^f^A^ + + 



+a?s|(^)=[I{iA 7 <»)'(»»A 7 (f<°' - aa 7 <»m°>} 

+i{iA 7 ^> - 2AH- M " + A 7 <?<<°>A a <f>>} 



+Aa<f < 2) - 2A,Ao5<" M1) + A 7 !r M0) A4f •'"] , (16) 

where S e is the spherical factor characteristic of iV-dimensional regularization. Here the hat 
indicates that all quantities are unrenormalized with respect to coupling constant and operator 
renormalization. The algebraic structure shown by the expression above follows from the renor- 
malization group. In addition to the anomalous dimensions one also encounters the coefficients 
of the beta-function. For instance (3q (see Eq. ([10])) is the lowest order coefficient, which also 
appears in the coupling constant renormalization, given by 

a s = a s (ii 2 ) 1 + a s (fi 2 )S e {2(3 ^-} , a s = ^ . (17) 

From the expression above one can in principle extract the first and second order anomalous 
dimension of the local operators in Eq. (|H]) which are given by A'y^'^ and A'y^'^ respectively. 
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However in the renormalization of the OME's one has to deal with two difficulties. The first 
one is caused by the fact that usually the external momentum p is taken off-shell (p 2 < 0). 
This means that the OME in Eq. fll5|) ceases to be a genuine S-matrix element and it becomes 
gauge dependent. Therefore one also has to carry out gauge parameter renormalization. The 
second problem, which is characteristic of spin operators, is the appearance of the 75-matrix 
and the Levi-Civita tensor e^ va p in the operator vertices (see Appendix A in U). In the case of 
iV-dimensional regularization one has to find a suitable prescription to define these essentially 
four dimensional quantities in iV-dimensions. In Jl| and || the vertex 7^75 is replaced by 

7^75 ^W^VV, (18) 

so that only the Levi-Civita tensor appears in the OME's. This prescription is equivalent to the 
one given by 't Hooft and Veltman |4| which is worked out in more detail by Breitenlohner and 
Maison || (HVBM). For the replacement of the 75-matrix in Eq. ([18]) see |6|], J7|. Although 
this prescription preserves the cyclicity of the traces it destroys the anticommutativity of the 75- 
matrix. This will mean that some Ward-identies or theorems will be violated. For example the 
non-singlet axial vector current, presented by the operator 0^ s '^ in Eq. (|14"D, gets renormalized 
in second order although it is conserved. Furthermore the Adler-Bardeen theorem |§ concerning 
the nonrenormalization of the Adler-anomaly is violated. This will affect the renormalization of 
the singlet axial- vector operator Ojp 1 ) in order a^. To undo these effects one has to introduce 
an additional renormalization constant in order to obtain the correct anomalous dimensions in 
the MS-scheme. The latter have now to be extracted from the renormalized rather than the 
unrenormalized OME's. After coupling constant renormalization the OME's are renormalized 
as follows 

INS,(n) = ~5,NS,(n)/~-l\NS,(n) JNS,(n) 

qq qq \ > qq qq 

A Un) = z ^ {z ^^n ) + {2 ^n)^^ (lg) 

where we have chosen the MS-scheme. In this scheme the constant for the HVBM-prescription 
can be written up to order a 2 as 



1. 



Z£.M(i) = 1 + a s S e [*&>.W] + a lSl [~A# (1) + <f ),(2) ] , (20) 



with r = NS, S. Notice that the difference between the singlet (S) and the non-singlet (NS) 
expression for shows up for the first time in second order (see |7|] ) . In this reference has 
been calculated for the first moment (n = 1) up to order a^. Recently this constant has been 
computed up to second order in || but now for general moments. In the non-singlet case it can 
be computed from the ratio 

7 5,NS,(n)/l\ _ A^g'( n \—p 2 ' / jj 2 , 1/g) | naive 

£ A qq ; (-p 2 //! 2 , 1/e) |hvbm 

where in the numerator one has used the so called naive prescription in which the 75-matrix 
anticommutes with all other 7-matrices irrespective of the value for the dimension N. The use 
of the naive method implies that the numerator can be replaced by the spin averaged OME 
A^'^ in which the 75-matrix does not appear. A similar derivation exists for Z^'^ where 
one also makes a comparison between the naive 75 and the HVBM prescription. From the 
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considerations presented above one could have asked the question why it is preferable to choose 
the HVBM instead of the naive prescription since in the latter case Z qq r '^ = 1? The reason is 

that the Levi-Civita tensor appears in the OME A gq which induces in the subgraphs containing 
quark lines the HVBM prescription. Therefore the naive method is inconsistent and it is better 
to use a consistent procedure like HVBM where all constants are fixed once and for all. The 
operator renormalization constants in Eq. (0), presented in the MS-scheme, read as follows 
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(22) 



The above expression differs from the usual one by the appearance of the term zjy with 
k = j = q which only contributes to {Z~ l ) qg (plus sign) and (Z~ l ) gq (minus sign) up to order 
ot 2 s . If this term is omitted then the anomalous dimensions will equal those present in Eq. (16|), 
which differ by a finite renormalization from the ones presented in the MS-scheme. Notice that 
the lowest order coefficients A^^'^ are not affected by any 75 prescription. Finally we want to 
emphasize that due to the pole term in Eq. fl2"0|) Z qq r does not represent a finite renormalization 
constant in the usual sense. Using the above procedure one can write the following expression 
for the renormalized OME in the MS-scheme 
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Notice that the coefficients Aa},- 



= («).(*) 



T6|). From the relation 



in this expression differ from the Aa^ present in Eq. 



A/fV)=4 n) (^")A/i 



f(«), 



(24) 



and Eq. ([TT] ) one concludes that the renormalized satisfy the renormalization group equation 
given by 



kj 



(25) 



After having discussed the renormalization of the OME's we now explain the procedure to 
compute the coefficient functions appearing in the spin structure function of Eq. (g) . They are 
obtained from the partonic subprocesses denoted by 
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+ i ->' X' . 



(26) 
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Figure 2: Contributions to the proces 7*+g — X' contributing to the partonic structure function 

9l,q- 



where i stands either for a quark (q) or a gluon (g) and 'X' represents an inclusive multi-partonic 
state. The above processes have been calculated up to order a 2 s in |J. Some Feynman graphs 
are depicted in Fig. ^ (i — q) and in Fig. 3 (i = g). The computation of the QCD corrections 
reveal ultraviolet, infrared and collinear divergences which appear in the Feynman integrals and 
in the phase space integrals. Like in the case of the operator matrix elements we regularize 
them by means of iV-dimensional regularization. Further we use the HVBM prescription for the 
7 5 -matrix as discussed below Eq. (|TS|). Adding all contributions one observes that the infrared 
singularities cancel and the radiative corrections are described by the parton structure functions 
which in general can be presented by the expression 
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Figure 3: Contributions to the proces ■y*+g — >' X' contributing to the partonic structure function 
9x,g- 
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with r = NS, S. The above expression still contains ultraviolet and collinear divergences both 
represented by the pole terms (l/e) k . The former are removed by coupling constant renor- 
malization (see Eq. fllTD). The residues of the collinear divergences are usually denoted by 
the DGLAP spin splitting functions AP^- which are related to the anomalous dimensions via a 
Mellin transform i.e. 



A 7; 



(2f 



The collinear divergences are removed by applying mass factorization which proceeds as follows 
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and the transition functions (r 



_l\r,(n) 
)ij 



are given by (see Eq. Q2"2])) 
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Using these expressions the longitudinal spin coefficient functions can be written as 
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so that they satisfy the renormalization group equations in Eq. dl2|) . Notice that in the above 



expression the coefficients Ac\( differ from Ac^ / in Eq. From the discussion above 

one infers that the renormalization of the operator matrix elements determine the way one 
has to perform the mass factorization on gl'^ and not vice versa. The reason is that the 
renormalization of the former (but not of the latter) is ruled by the Ward-identities and some 
theorems which are violated by the HVBM-prescription. This has forced us to introduce the 
additional constant Z^ r '^ in Eq. (^Op to restore the wanted properties on the level of the 
renormalized operator matrix elements presented in Eq. (|23|). If we would have accidentally put 
^qq = 1 the coefficient functions and the renormalized operator matrix elements get different 
anomalous dimensions and the relations in Eq. ( |30|) would be violated. From Eqs. fl24|), fl25|) it 
also follows that the parton densities and the coefficient functions would have different anomalous 
dimensions. Hence the structure function gi(x,Q 2 ) @ would not satisfy Eq. ( |T3"D anymore so 
that it is no longer a physical quantity. Therefore this choice for is unacceptable. 
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2 Fragmentation into polarized hadrons in electron- 
positron annihilation 

Single hadron (denoted by H) inclusive production in electron-positron annihilation is given by 
the reaction 

e-(k, ^) + e+(Z 2 , a 2 ) - V(q) - H(p, s) +' X' . (32) 

Here we have introduced a similar notation to the one in reaction ([!]). However the Bj0rken 
scaling variable is defined by 

x = , q 2 = Q 2 >0, < x < 1 , (33) 

Q 

for timelike momenta of the vector boson V. The annihilation process can be depicted as in 
Fig. [I] where now the incoming hadron is outgoing and the electron in the final state becomes a 
positron in the initial state. In the case the incoming electron in reaction (|32| ) is longitudinally 
polarized downwards, i.e. o\ =1, and the positron is unpolarized then one can simplify the cross 
section when the process becomes purely electromagnetic. In this case V = 7 and we get 

da n ^(l) da u ^(l) _ N 71a 2 ^ H q2 
dxdcos8 dxdcos8 Q 2 ' 

The above expression represents the difference between the cross sections where the detected 
hadron H is polarized parallel s =| or anti-parallel s =| with respect to the electron. Further 
Nc denotes the colour factor in SU(Nc) and 6 is the polar angle describing the direction of 
the momentum of H in the CM. frame with respect to the incoming electron. Notice that the 
hadron fragmentation function gf(x,Q 2 ) can be also measured in unpolarized electron-positron 
scattering provided reaction ([$2]) is mediated by the Z-boson. Here it appears via the axial- vector 
coupling of this boson to the lepton-pair. The above hadron fragmentation function shows up 
in the anti-symmetric part of the hadronic structure tensor 

W^(p,q,s) = -Ly d*ze**(0 I J,{z) \ H(p, s), X)(H(p, s), X | J„(0) | 0) , (35) 

which can be decomposed in the same way as shown for Eq. @ in deep-inelastic scattering. 
Since H is exclusive the above expression is not a Fourier transform of a curent-current correlation 
function which implies that we cannot insert the lightcone expansion of Eq. (|15[) . Therefore we 
can only work in the QCD improved parton model picture in which the fragmentation function 
has the following form 



Uf Jx z 1 H K z ' ' yv n z ' 

+ADM(|, „»)<,(*, J) + „ y A D r S (^ 2 )AC(^f )] ■ P6) 

The spin parton fragmentation densities denoted by ADf^(z, fi 2 ) are the analogues of the parton 
densities in Eq. (||) and they satisfy the same renormalization group equations. However beyond 
lowest order the anomalous dimensions are different for these two densities (see [ID], [|TI|, [T2||). 



The anomalous dimensions ruling the evolution of the spin fragmentation densities have recently 



S 



been calculated up to second order in fl2 |. To obtain the timelike spin coefficient functions one 



has to calculate the timelike photon analogues of the graphs in Figs. |^, 3. Here the incoming 
quark and gluon become outgoing and they now fragment into the hadron H. Furthermore the 
spacelike photon turns into a timelike one. The calculation of these coefficient functions was 



recently done up to second order in [jn|. It proceeds in the same way as for deep inelastic 



scattering in section 1 where again the HVBM prescription for the 75-matrix is chosen. After 
having calculated the parton fragmentation functions denoted by gf'f one has to perform mass 
factorization analogously to Eq. (fffi). However our calculation reveals that the renormalization 
constant Z qq r is different for timelike (fragmentation function) and spacelike (structure function) 
quantities. In one has obtained the non-singlet part of this constant for the timelike (T) 
process fl3"2|) by computing the ratio 

7 5,NS,T,(n) _ yi£ W /P ; L / b ) I /o 7 \ 

^NS,(n) (g2/ ^ 1/£) y=Q>> 

where ^^ NS, ( n ) j s the parton fragmentation function in unpolarized electron-positron annihila- 
tion. It arises due to the axial-vector coupling of the Z-boson to the outgoing quark anti-quark 
pair (see Fig. |l|). If the 75 anticommutes with the other 7- matrices one can show that the 
coefficient functions ACf^ and C^'^ are equal up to order a 2 . Since the HVBM prescription 

destroys this property for gfg NS '^ n ' ) we have to multiply the latter by Z qq NS,T in order to obtain 
the correct coefficient functions. If we assign to the constant in Eq. (|20"D the superscript S (here 
spacelike), the following relation holds for the inverse Mellin transforms 

Z 5 qq NS ' T (z) = -zZ^i 1 -) + a^^Hz) In,] . (38) 



The above equalitity demonstrates the breakdown of the Gribov-Lipatov relation J14 in order 
a 2 s . The above relation is also found for the timelike and spacelike non-singlet splitting functions 



P qq in unpolarized scattering in Eqs. (6.37) and (6.38) of |L0| where z qq > is replaced by the 
lowest order DGLAP splitting function P^'^. Notice that the same relation holds for the spin 
splitting functions because AP^ = P^ (see |12|]). The dependence of Z qq NS on the quantity 
under consideration reveals that it is not an universal constant. Therefore aside from the pole 
term already mentioned above Eq. (^), it does not represent a genuine renormalization constant 
in the usual sense. 
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